Trial wavefunctions for the Goldstone mode in v — 1/2 + 1/2 quantum Hall bilayers 
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Based on the known physics of the excitonic superfiuid or 111 state of the quantum Hall v — 
1/2 + 1/2 bilayer, we create a simple trial wavefunction ansatz for constructing a low energy branch 
of (Goldstone) excitations by taking the overall ground state and boosting one layer with respect to 
the other. This ansatz works extremely well for any interlayer spacing. For small d this is simply the 
physics of the Goldstone mode, whereas for large d this is a reflection of composite fermion physics. 
We find hints that certain aspects of composite fermion physics persist to low d whereas certain 
aspects of Goldstone mode physics persist to high d. Using these results we show nonmonotonic 
behavior of the Goldstone mode velocity as a function of d. 



The v = 1/2 + 1/2 quantum Hall bilayer is a remark- 
ably rich system [2 H] . At small enough spacing between 
the layers, d, the system is known to be an excitonic su- 
perfiuid [3] known sometimes as the 111 phase [I]. At 
larger layer spacing, a phase transition or crossover is 
observed experimentally [51-11 1 j leading to a compress- 
ible phase which is well described by two weakly cou- 
pled composite fermion Fermi liquids. The nature of this 
crossover, and whether there are intervening phases be- 
tween small and large d, has been a matter of some debate 
in the community [T2TU9j . 

There are some results, however, that are extremely 
well established theoretically. In the limit where d be- 
comes very small, it is known that the Halperin 111 trial 
wavefunction becomes exact pE]. In a more BCS-like lan- 
guage this wavefunction can be expressed as [3] 

ini)=n(4t+ c L)io> a) 

X 

where t an d I indicate the layer index (we assume the 
real spin is frozen throughout this paper) and X con- 
stitutes the orbital index within the Lowest Landau level 
(chosen to be the ^-directed momentum in Landau gauge 
for example) |20] . 

The BCS-like form of Eq. [I] allows one to consider long 
wavelength Goldstone excitations of the form [3l [21] 

|111 - excitation, k) = ]J( c x,T + e%kX c x,i)\ Q ) ( 2 ) 
x 

These modes are expected to form a linearly dispersing 
low energy branch with energy proportional to k for small 
k. Physically, this Goldstone mode corresponds to super- 
flow — one layer being boosted with respect to the other. 
Both a linearly dispersing mode }22] and excitonic super- 
flow [71 [SI |2"3T - T2"5] were observed experimentally in this 
system. 

Away from the d — > limit the form of the bilayer 
ground state is not known exactly. However, so long as we 
remain in the same phase of matter, there will continue 
to be a linearly dispersing Goldstone mode in the long 



wavelength limit. An approximate expression for this 
Goldstone mode can be obtained from the ground state 
wavefunction at any d simply by boosting one layer with 
respect to the other. One purpose of the current paper 
is to test this technique of generating trial wavefunctions 
for the long wavelength Goldstone modes. 

We note that this technique is not expected to be exact 
away from d = 0, but for small d is expected to be quite 
accurate. In a conventional picture of superfluidity, one 
might imagine that it would be better to find a way to 
boost the superfiuid fraction while leaving the "normal" 
fraction unboosted (only at d = is the system entirely 
super in some sense |13j). Nonetheless, our technique 
appears to work quite well even away from d = 0. 

One might expect that once the system is no longer 
in the 111 phase of matter (roughly d > 1.5 mag- 
netic lengths), our technique for generating excited states 
would fail. However, this turns out not to be the case. 
First of all, at intermediate d there may exist an inter- 
layer paired state as discussed in Ref. [TBI Such a paired 
state would also have a Goldstone mode that could be 
generated from the ground state by boosting one layer in 
exactly the same way. 

However, even at very large d when such a pairing 
phase is either absent or pairing is extremely weak, our 
scheme for generating excited states still works surpris- 
ingly well. To understand why this is so, we realize that 
at large enough d each layer is essentially independent. 
To first approximation, each layer forms a composite 
fermion Fermi liquid, which for finite size system has fi- 
nite momentum (except when the number of electrons 
exactly fills a shell). The two Fermi liquids are weakly 
coupled and can combine their momenta to form an over- 
all zero momentum ground state. But since the coupling 
between the two layers is weak, it costs very little en- 
ergy to instead form a state of overall finite momentum 
— which can be interpreted as boosting one layer with 
respect to the other in comparison to the ground state. 

In the absence of any superfiuid order parameter (at 
large d) it is probably not strictly appropriate to refer 
to this low energy mode as a Goldstone mode. However, 
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FIG. 1. Overlaps of the trial states |3| with the exact eigenstates of the Coulomb Hamiltonian for a selection of system sizes 
N = 10, 12, 14 and 16 on the sphere. The magnitude of the overlap of the trial state with the exact groundstate is indicated by 
the size of red dark shaded circles. Blue lightly shaded circles additionally indicate the overlap of the same trial states with the 
first excited state of the exact spectrum. Overlaps at L = are equal to one by definition and give the overall scale. The trial 
states are very accurate at d = 0, giving a good description of the lowest energy mode in each sector of angular momentum 
(the Goldstone mode) up to large L. At finite layer separation d, the Goldstone mode is always present at small L, but does 
not reach to similarly high values of angular momentum. The description is again more successful at very large d. 



since this mode may evolve continuously into the Gold- 
stone mode at smaller d we will abuse nomenclature and 
continue to call it a Goldstone mode. (If, as conjectured 
in Ref. HH the bilayer is actually paired out to large d, 
then the usage remains correct). 

Throughout this paper we will work with a spherical 
geometry. In this case, boosting one layer with respect 
to the other corresponds to applying the angular mo- 
mentum raising operator L + to one layer but not the 
other (call this operator meaning that it is applied 
to the f layer only). In the appendix we show that if we 
start with any L = state of the entire system, apply- 
ing (L + ^) J generates a bilayer state with overall angular 
momentum L = L z = J . Our technique is then to use 
exact diagonalization to generate the L = ground state 
of the bilayer system, which is used to obtain the trial 
wavefunction for the excited state 

| Trial(d) : L = J) = (L+ T ) J | Ground State(d) : L = 0). 

(3) 

In turn, we compare this trial state to the exact excited 
states with angular momentum L = J. 

Our numerical work is based on exact diagonalization 



of the Coulomb Hamiltonian for a bilayer system on the 
sphere [35]. We simplify the problem to exclude issues 
related to a finite tunneling amplitude between the lay- 
ers, Landau-level mixing, or spin (which we assume is 
polarized), and model each layer as an ideal 2D plane 
without considering its width into the third dimension. 
At fillings smaller than one per layer, the Hamiltonian is 
thus given by the projection of the Coulomb interaction 
into the lowest Landau-level 



(4) 

where sums run over all particles with the given pseudo- 
spin. The interactions are parametrized by the layer sep- 
aration d that is measured in units of the magnetic length 
£q = y/hc/eB. All lengths given in this paper should be 
understood to be measured in units of Hq, where this is 
not explicitly indicated. 

In our exact diagonalization calculations, we focus 
specifically on the density-balanced bilayer system with 
iVf = iVj. = N/2 and consider the state at the shift of 
the 111 state, namely N$ = N — 1. We obtain the two 
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FIG. 2. Overlaps of the trial states for the Goldstone 
mode at d — (top) and d — 3 (bottom) with the exact 
state, as in Fig. [Tj plotted here as a function of wavevec- 
tor k ~ L I \J~N~^fl. Despite the size of the Hilbert spaces 
increasing significantly between the smallest and largest sys- 
tem shown, the overlap remains roughly constant or maybe 
slightly increases with N. The region of high overlap extends 
roughly up to k w 2i^ x . The failure of this approach for d = 3 
and N = 12 is discussed in detail in the text. 



lowest-lying eigenvalues and eigenvectors in each sector 
of angular momentum. This is most easily achieved using 
a projected Lanczos algorithm |27] which uses an addi- 
tional projection to the lowest energy subspace of min- 
imal angular momentum after each multiplication with 
the Hamiltonian. In a given sector with fixed L z , this 
procedure therefore directly yields eigenstates of L 2 with 
the eigenvalue L = L z . 

To assess the accuracy of the trial states (|3| for the 
Goldstone mode, we consider their overlap with the low- 
est energy state in each sector of angular momentum L 2 
and for layer separations d = . . . 3£q, in steps of \l§ [28] , 
We note again that the the trial states are generated by 
applying the operator (L + ^) J to the exact groundstate 
at L = 0. Only at d = is an exact analytical expression 
of the groundstate known - the Ill-state. At other values 
of d, the numerical groundstates from exact diagonaliza- 
tion are used, although very accurate trial wavefunctions 
are also known [19j . The results are summarized in Fig.JTj 
which also indicates overlaps with the first excited state 
in addition to the overlaps with lowest energy state in 
each sector. At d = 0, the ansatz ^ is very successful, 
describing excited states up to high angular momentum 
accurately. This is shown in more detail in Fig. [2] a), 
which displays the magnitude of the overlap as a func- 
tion of wavevector k. Surprisingly, for d = the overlap is 
very consistent with system size at given k, even though 
the Hilbert space dimension increases strongly with TV. 
A very good description with overlaps above 0.8 is given 
up to a wavevector of k ~ 2£q 1 . Turning back to Fig. [I] 
we now focus on the overlaps at finite values of layer sep- 



aration. At d — 0.5£q, our trial states obtain significant 
overlaps only with the first excited state at L = 1. This 
is not due to a disappearance of the linearly dispersing 
mode however. Rather, a level crossing appears with dis- 
tinct excitations occurring at energies less than that of 
the Goldstone mode, as can be seen from the significant 
overlaps with the first excited excited state in the sector 
of L = 2. At d > 1.54) this overlap with the first excited 
state again disappears, signalling the presence of addi- 
tional low-lying excitations of a nature different from the 
linearly dispersing Goldstone mode. 

Finally at the largest value of d = 3£q, the ansatz for 
the boosted trial wavefunctions becomes more accurate 
than at intermediate d, signalling the possible emergence 
of a distinct mode of low-lying excitations. Based on the 
overlap with the trial states, we can point out that there 
are strong finite size effects in the physics at large d. As 
with d = 0, Fig. [2] b) displays the numerical values of 
the overlaps at d = 3. These data single out the system 
with N = 12 particles as particularly poorly described 
by these trial states. Here, the physics at large d is 
clearly dominated by the shell filling effects of composite 
fermions. As N-^ — — 6 electrons per layer precisely 
fill the lowest two shells of composite fermion orbitals in 
one quantum of effective flux, this system size is aliased 
with the situation where each layer forms its own incom- 
pressible v = 2/5 state with angular momentum zero in 
each layer — making it impossible to form higher angu- 
lar momenta states by boosting one layer with respect to 
each other. 

Before proceeding further, note the rather unusual fea- 
ture that, excepting N — 12, the trial states give higher 
overlaps for larger systems at d = 3£q. This unusual 
behaviour is related to a different manifestation of the 
composite fermion shell filling effect. As we have shown 
in Ref. [TS1 the groundstate at large layer separation is 
a state in which each layer individually obeys Hund's 
rule and maximizes the angular momentum per layer 
L a = =L\, while both layers are combined into a 
total L = state. Without modifying the correlations 
inside each layer, the same states with L a per layer can 
be paired into excited states with subsequently larger an- 
gular momenta, up to a maximum L max = 2L a . For 
the system sizes with partially filled composite fermion 
shells in each layer, one obtains the values of L a = | for 
N = 10, L a = | for N = 14, and L a = § + § = 4 for 
N = 16 particles. We therefore expect a low-lying mode 
of excitations with angular momenta up to L max = 3, 
L max — 5, and L max = 8, respectively. Indeed, upon 
inspection of the spectra, such a mode can be identified. 

As an example, Fig. [3] displays the spectra for the 
system with N — 16 particles, including different val- 
ues of the layer separation. Indeed for d — 3£q, shown 
in the bottom right panel, there is clear evidence of a 
mode of excitations terminating at L = 8. Its disper- 
sion is approximately linear at small L, however, it has a 
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FIG. 3. Spectra of the Coulomb Hamiltonian in the bilayer system on the sphere for N = 16 particles with N$ = 15 flux, for 
layer separations d ranging from d = (top left) to d — 3£o (bottom right). Energies are indicated in units of e 2 /e£o relative 
to the groundstate energy _Egs at L = 0. Red crosses mark the lowest two eigenstates in each sector of given L. Blue circles 
indicate those states which were identified as part of the Goldstone branch by high overlaps (see Fig. [I]). At d = the Goldstone 
mode is very well formed, while it is clearly visible at small and intermediate d how there are level crossings with additional 
low lying excitations. At d — 3, a mode of low-lying excitations is once again clearly separated from the rest of the spectrum, 
which terminates sharply at L = 8. Its dispersion is linear at small L, but it has a quadratic component responsible of the 
upturn at larger L. Note the change in the scale of the y-axis for the three columns of panels. 



quadratic component as well. This should be compared 
to the Goldstone mode at d = (top left panel), for which 
linear dispersion is clearly realized up to high values of 
angular momentum. 

Once the termination of the low energy branch at 
L max = 8 is identified at large d, it becomes apparent 
that this feature of a jump in the spectrum at L max exists 
at all values of layer separation shown, with the exception 
of the SU(2) invariant case of vanishing d. Note that this 
termination is a feature of composite fermion physics, 
explained by successively filling the lowest shells of these 
composite particles, while obeying Hund's rule. The ob- 
servation that composite fermion physics intervenes at 
very small layer separation had been made previously by 
the current authors. While the Ill-state can be regarded 
as a condensate of composite bosons, it was shown that an 
accurate description of the groundstate requires a mixed- 
fluid description of both composite bosons and composite 
fermions at any finite layer separation [IB"], H~9] . The iden- 
tified jump may constitute evidence for the mixed-fluid 
picture in the excitation spectrum, but more study of 
these excitation will certainly be required. 

A linearly dispersing mode at small k exists at all val- 
ues of the layer separation. The states which were shown 
in Fig. [I] to have large overlap with the trial states ([3| 
are highlighted by blue circles in the spectra shown in 
Fig. [3j These states very accurately come to lie on a 
single line, which is true especially for the Ill-state at 
d — 0, but also for the intermediate layer separations 
such as d = 1£q, where the first excited state at L = 2 



lies in the continuation of the line through the points at 
L = and L = 1 and is shown to be associated to the 
Goldstone mode by its overlap. Judging by the spectra, it 
is also very suggestive that multiple level crossings occur 
at larger values of L, for example at L = 2 for d = 1.5£o- 
Finally, between d = 2£q and d = 3£q, a change occurs 
in the association of the low-lying mode with the Gold- 
stone mode trial states. For d = 2£q only the states up 
to L = 2 have a good overlap, and the remaining states 
of the already well formed band of low-lying states in the 
exact spectrum are of a different nature. The transition 
to low overlaps occurs at a point where this band has a 
visible kink and flattens out. Finally at d — 3£o, this 
low lying band has good overlaps up to much higher mo- 
menta, which, as we have discussed above, is a reflection 
of two approximately uncoupled composite fermion Fermi 
seas which each maximizes its own angular momentum 
according to Hund's rule. 

Given the existence of a linearly dispersing mode over 
the whole range of layer separations, we now consider 
how its velocity changes with d. Comparing the results 
obtained for different system sizes, a relatively strong 
dependence of the velocity v = dE/dk is evident. Ap- 
plying charging corrections to take account of the shift 
in the charge-density of the system [29 does not suf- 
fice to absorb these effects. Thus, in addition to mea- 
suring energies in units of the rescaled magnetic length 
£' = y/vN^/NtQ, we analyze the scaling of the mode 
velocity as a function of the inverse system size N —1 . 
These scalings are fitted well by linear extrapolation to 
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FIG. 4. The velocity of the Goldstone mode depends strongly on the layer separation, and peaks near d ~ 1$. The velocity 
has significant finite-size effects on the sphere, therefore we extract an estimate of its value in the thermodynamic limit by 
extrapolating the finite size values over the inverse system size (see inset). 



the thermodynamic limit, as shown in the inset of Fig. [4] 
The data for N = 12 and large layer separation are eas- 
ily identified as outliers, due to the shell filling effects 
discussed above. Generally, the slope is decreasing with 
system size. However, at small finite layer separations 
the opposite scaling takes effect. The resulting estimates 
for the mode velocity in the thermodynamic limit shown 
in the main graph of Fig. [4] therefore show a pronounced 
maximum near d = £q, which has about twice the mag- 
nitude as the value found at d = 0. Beyond this point, 
the mode velocity decreases monotonically and drops to 
about l/10 th the value of the Ill-state at layer separa- 
tion d = 3£q. 

A previous experiment [55] has probed the velocity of 
the neutral mode at layer separations near the transi- 
tion into the incompressible phase at large d, obtaining 
a combined best fit of v = 1.4 x 10 4 m/s for data at 
the three layer separations d\ = 1.614), <^2 = 1-714 
and c?3 = 1.764- Based on linear extrapolation be- 
tween our numerical data at d = 1.54 and d = 24, 
the corresponding estimates are v(di) — 1.14 x 10 4 m/s, 
v(d 2 ) = T21 x 10 4 m/s, and v(d 3 ) = 1.33 x 10 4 m/s, all 
slightly smaller but within about 20% of the proposed 
fit to the experimentally obtained values. Had the data 
in Ref. 22]been fitted separately at each layer separation, 
the velocity at d = 1.714 would be estimated to be about 
10% smaller than that at d = 1.614), roughly reflecting 



the ratio of our predicted values. The data at d — 1.764 
appears to be rather noisier, probably due to the vicinity 
to the phase transition, and would be difficult to fit on 
its own. We suggest that a significant enhancement of 
the linear mode velocity should be seen deeper inside the 
inter-layer coherent phase at smaller layer separation. 

To summarize our results, we use the ansatz Eq. [3] 
to construct trial wavefunctions for a low energy branch 
of excitations based on the exact ground state wavcfunc- 
tion. This ansatz is accurate at all interlayer spacings d 
when k is small, and it is accurate at all k when either d is 
small or d is large (so long as we do not have a filled shell 
configuration, whereupon only k = is in this low energy 
branch). We find hints that certain aspects of the com- 
posite fermion physics persist to low d whereas certain 
aspects of the Goldstone mode physics persist to high d. 
Applying these results to the analysis of our numerical 
data we show nonmonotonic behavior of the Goldstone 
mode velocity as a function of the layer separation d. It 
would be interesting to look for this nonmonotonicity of 
the Goldstone mode velocity experimentally. 
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bridge, as well as from the Institute of Complex Adap- 
tive Matter (ICAM). SHS thanks the Aspen Center for 
Physics for its hospitality. 
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Appendix: Angular Momentum of Two Coupled 
Subsystems 

In this appendix, we will use the standard angular mo- 
mentum notation \l, m) to indicate eigenstates of the L 2 
and L z operators. 

We consider two subsystems f and 1 with correspond- 
ing angular momentum operators Lq and Lq. These two 
subsystems combine to form the total system with angu- 
lar momentum operator 



J 4 



(5) 



Our objective is to show that given an eigenstate of the 
total system with \l = 0, m = 0) application of (i+f) J to 
this system will produce an eigenstate of the total system 
with \l = J, m = J). To achieve this, it is sufficient to 
show that 

L +t \J,J)~\J+l,J+l) 

Obviously applying to any state increments its over- 
all L z eigenvalue by one (or kills the state), so all that 
remains is to show that applying to | J, J) results in 
an eigenstate L = J+ 1 of L 2 , i.e, results in an eigenvalue 
of L 2 being given by (J + 1)( J + 2). 

Using Eq. [5] it is just a matter of some algebra to show 
that 

[L+f, L 2 ] = 2L z ^L + i - 2L + ^L z _i) 

= 2(L+-L +t )(L z + l)-2L zl L + 

We then apply both sides of this equation to the state 
| J, J) . Noting that L + kills | J, J) , we obtain 

L 2 [L +t \J,J)} = [L +t (J(J+l)) + 2L +t (J + l)]\J,J) 
= (J + l)(J + 2)[L +t |J,J)] 

which completes the proof. 
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